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this difficulty can be lessened ; by taking 500 days instead of iooo 
as the second unit. But this will detract from simplicity of form* 

A yet more serious drawback to the theoretical investigator is 
that the fundamental epochs usually adopted in astronomy* and 
for which the elements must be found, do not correspond to- any 
power of 10 in the days of the Julian period. A complete trans¬ 
formation of the elements is therefore required to form the numbers 
on which the tables are based. If, therefore, multiples of 500 or 
1000 days are used instead of years, I should prefer to count them 
back from 1900*0, thus gaining all the advantages of the Julian 
period without any other disadvantage than that of non-correspond¬ 
ence with the eclipse and other tables of Oppolzer. 

The reduction of such a system to the ordinary calendar may 
be made a very simple matter. It seems to me, therefore, that the 
maximum of advantage will be reached by giving the fundamental 
arguments for cycles and periods based on multiples of 500 days 
before and after the fundamental epoch 1900 Jan. o. 

Probably the most convenient fundamental quantities to tabu¬ 
late will be the longitude of the node, and the mean distances of 
the Moon and of its perigee from the node, all expressed in circum¬ 
ferential units. Then, whatever form the tables may be thrown 
into, we shall have the nearest approach to a simple, straight¬ 
ahead computation. 

Finally, a serious problem is that of summing perhaps 100 
periodic terms with coefficients not differing greatly from o"*oi. 
I have devised a machine for this purpose, the description of which 
must form the subject of another publication. 


An Example of Professor Karl Pearson’s Calculation of Correla¬ 
tion in the case of the Periodic Inequalities of Long-period 
Variables. By H. H. Turner, D.Sc., F.R.S., Savilian 
Professor. 

1. The following note is written with a twofold purpose, 
firstly, it is hoped that an indication of some value has been 
obtained with regard to the features of “long-period” variability; 
and secondly, the opportunity is taken to write out in full a simple 
example of the calculation of “ correlation ” between quantities by 
the methods of Professor Karl Pearson. 

In the M.N. for March last (p. 416) Professor Pearson himself 
gave au admirable summary of methods; but he naturally did not 
repeat the elementary working which has become so familiar to 
him, and has been given often before in other connections. There 
are doubtless mauy to whom this working is already familiar; 
but there are certainly many others who do not know it and who 
might use it if they had an astronomical example readily accessible. 
In these busy days many people have not the leisure to search 
for references in scientific literature outside their own subject. 
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r 2. In support of pay* view that there is a'esed; fc# amexample, 
I may cite an illustration shown me by Professor;!*}. ,C. Pickering 
(see Observatory , March 1905, p. 153). 

Suppose we have the following observed values of A and B 


A 

B 


5 

6 


6 

6 


7 

6 


6 

7 


7 

7 


8 

7 


7 

8 


8 

8 


9 

'8 


v First group according to B, in sets of three for which B has 
the same value. 


Then when 


Therefore 


B = 6 7 8 

Mean A = 6 7 8 

A = B. 


Secondly, group according to A. 

Then when A=5 6 7 8 9 

Mean B = 6 6 \ 7 7J 8 

Therefore (A - 7) = 2(B - 7). 

How, neither Professor Pickering nor myself was prepared at 
that time to deal with this situation by a definite process; and 
yet this is just an elementary case of the kind which Professor 
Pearson’s methods were devised to meet. 

This example is worked out below, and it is shown that the 
numerical measure of the correlation is , 

r=*o8 + *22 

or almost nothing at all: so that, in spite of appearances to the 
contrary, we are not entitled to assume any relationship between 
A and B. Putting it in another way, one proposed relation 
(A = B) is as good as another (A — 7) = 2 (B — 7). 

3. A point of detail may be mentioned here. In much 
statistical work, a large number of figures are used. Thus we get 
such statements* as 


r = 


1262*51 


2870 x 1*13637 x 1*52135 


*? 5445 - 


The probable error of r being about ±’01, some of these figures 
are superfluous; and in what follows fewer figures are used. 
This, however, represents a personal view which is, I find, not 
generally approved by other workers. 

4. The particular example selected for treatment is the 
discussion of the elements of maximum given by Chandler for 
long-period variables. The following particulars are taken from 


* Frequency Curves and Correlation , by W. P. Elderton, p. 119: an 
excellent little book, from which much is to be learned. 
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his “ Revision of Elements of Third Catalogue ” (Astronomical 
Journal , No. 553) :— 

Elements of Maximum. 

Star. Epoch. Period. Inequalities. 

R Androm. 2400141 + 410*3 E + 30 sin (12 0 E +90°) 

S Cassiop. 2401603 -b 610*5 E + 37 sin (15 0 E + 59 0 ) 

and so on. The mean period of S Cassiop. is 610*5 days; but 
there is an inequality which displaces the maximum 37 days one 
way, and 12 periods later (since 180715° = 12) the displacement 
is 37 days in the other direction. At these times the period is 
about the mean, but midway between such times it is longer and 
shorter. In the “Revision” above mentioned, Chandler gives 
such inequalities for 37 stars, as shown in Table I. We shall 
for the present only consider the period P, the coefficient C, and 
the argument A; thus for S Cassiop. P= 610*5, C = 37, and 

A= 15 0 . The epoch 59 0 does not concern us at present. 


Table I, 


Star’s Name. 

No. 

P 

c 

A 

Phase. 

M - m 

aP 

R Androm. 

112 

410 

30 

12 

218 

120 

- 170 

S Cassiop. 

432 

610 

37 

15 

148 

280 

- 50 

S Piscium 

434 

404 

18 

IO 

. 57 . 

163 

78 

R Arietis 

782 

187 

7 

5 

292 

92 

- 3 

o Ceti 

806 

332 

18 

5 

... 

125 

- 82 

R Persei 

1222 

210 

15 

8. 

224 

96 

- 18 

R Aurigse 

1855 

459 

19 

12 

3 2 ° 

235 

+ 11 

R Lyncis 

2478 

379 

14 

15 

249 

186 

- 7 

R Geniin. 

25^8 

'370 

35 

6 

121 

121 

- 128 

S Can. Min. 

2684 

'330 

20 

12 

120 

164 

- 2 

R Cancri 

2946 

362 

60 

6 

34 i 

125 

- 112 

R Carinse 

. 34 i 8 

3 IQ 

25 

9 

289 

136 

- 38 

R Leo Min. 

3477 

37 i 

20 

10 

11 

165 

- 41 

R Urs. Maj. 

3825 

302 

1 r 

8 

63 . 

I IO 

- 82 

T Urs. Maj. 

: 45 H 

257 

20 

9 

203 

108 

- 42 

R Yirginis 

. 4521 

145 

20 

2 

103 

68 

- 9 

S Urs. Maj. 

4557 

226 

35 , 

5 

3°9 

108 

- 10 

S Bootis 

5157 

270 ■ 

60 

4 

■ 87 

132 

- 6 

R Camelop. 

5!90 

270 

65 

4 

244 

142 

+ 14 

R Bootis 

5237 

223 

9 

9 

251 

102 

- 20 
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Table I. — continued. 


Star’s Name. 

No. 

p 

c 

A 

Phase. 

M -m 

aP 

S Serpen tis 

55 °I 

369 

116 

0 

4 

0 

122 

... 

... 

S Coronae 

5504 

361 

8 

12 

8l 

120 

- 121 

R Serpeptis 

5677 

357 

35 

4 

IO9 

151 

- 55 

R Herculis 

577 ° 

318 

17 

10 

29 

... 

• •• . 

W Herculis 

595 ° 

280 

26 

13 

288 

128 

- 24 

S Herculis 

6044 

3°8 

35 

9 

239 

152 

- 4 

T Herculis 

6512 

165 

10 

5 

152 

79 

- 7 

R Sagittarii 

6905 

269 

18 

10 

18 

138 

+ 7 

S Sagittarii 

6921 

231 

i 5 

10 

167 

102 

- 27 

R Delphini 

726l 

284 ■ 

26 

9 

266 

130 

- 24 

U Capric 

7455 

203 

20 

5 

67 

... 

... 

T Aquarii 

7468 

203 

8 

7 

357 

88 

- 27 

R Yulpec. 

7560’ 

137 

18 

4 

127 

62 

- 13 

X Capric. 

7577 

218 

20 

10 

98 

11 7 

+ 16 

R Pegasi 

8290 

378 

60 

8 

7 i 

172 

- 34 

R Aquarii 

8512 

387 

35 

10 

97 

... 

... 

R Cassiop. 

8600 

432 

32 

9 

233 

182 

- 68 


5. The first two columns of Table I. give the star's name and 
its number on Chandler's system. The third column, P, gives the 
period in days; C is the coefficient of the periodic inequality in 
days, and A the coefficient of E within the bracket. The explana¬ 
tion of the last three columns will be given later. 

6. Now the example selected to be given in full is that of the 
relation between P and A. To determine whether there is such 
a relation, we could, of course, proceed in the ordinary way to 
group together stars with nearly the same P, and take the mean 
values of A for them. This would give us a fair indication, but no 
notion of its probable error. Professor Pearson's procedure is not 
much longer; it includes the ordinary procedure, and it gives us 
a definite measure of the probability of the relation being real. 

7. We first form a “correlation table” as in Table II. 


[Table v 
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Table IL 


Correlation Table for P and A. 


Pin days. \^A= 

2° 

,, 3° 

4 ° 

5° 

6° 

7° , 

8° 

9° 

IO° 

ii ° 

12° 

*3* 

1 4° 

15° 

Sums 

601 to 650 














I 

36 

I 

551 to 600 














O 

30 

0 

501 to 550 














O 

24 

O 

451 to 500 








O 

0 1 

_ 3 ' 

0 

6 

I 

9 

0 

O 

O 

18 

I 

401 to 450 








I 

I 1 

_ 2 

O 

4 

I 

6 

O 

_ 8 

i O 
_10 

O 

_ 12 

D 

351 to 400 



2 

O 

4 

2 

3 

O 

_? 

I 

1 

O 

2 

1 

O 

1 _ Q 

I 

O 

4 

1 O 

5 

I 

_ . _6 

9 

301 to 350 



O 

I 

0 

0 

I 

2 

I 

0 

I 

O 

O 

0 

6 

251 to 300 



2 

5 

O 

_ 4 

O 

3 

0 

_ 2 

O 

1 

2 

I 

1 

1 0 

2 

O 

I 

4 

1 O 

5 

O 

6 

6 

201 tO 250 



O 

10 

2 

8 

O 

6 

I 

4 

I 

2 

I 

2 

2 

1 O 

4 

O 

6 

O 

8 

1 0 

- 10 

O 

12 

7 

151 to 200 



O 

IS 

2 

_ 12 

O 

9 

O 

6 

O 

3 








2 

101 to 150 

I 

28 

O 

_2* 

0 

I 

_2P 

5 

O 

O 

12 

O 

_8 

O 








2 

Sams 

I 

! 5 

2 

I 

3 

6 

7 ‘ 

! 0 

4 

.1 . 

0 

2 

37 





v ,= 

|/ 2 = 




vi - 

i/ 2 = 



A 

Sum. 

Factor. Product. 

Product. 

P Sum. Factor. 

Product. 

Product. 



O 

2 

I 

-7 

- 7 

49 

125 

2 

“4 

- 8 

32 



3 

O 

-6 

# 

0 

0 

175 

2 

“3 

- 6 

18 



4 

5 

-5 

-25 

125 

225 

7 

- 2 

-14 

■28 



5 

5 

-4 

- 20 

80 

2 75 

6 

- 1 

- 6 

6 



6 

2 

- 3 

- 6 

18 

325 

6 

0 

0 

0 

Products S'(AP 

7 

1 

-2 

- 2 

4 

375 

9 

4 I 

+ 9 

9 

7 17 

+64 

8 

3 

- 1 

- 3 

3 

425 

3 

4-2 

4 - 6 

12 

+104 

~ 9 

9 

6 

0 

0 

0 

475 

1 

+ 3 ' 

+ 3 

9 


= + 142 

10 

7 

+ 1 

+ 7 

7 

525 

0 

+ 4 

0 

0 

37(--84H- 

•27)=+ 9 

11 

0 

4-2 

0 

0 

.575 

0 

+ 5 

0 

0 

Difference = 133 

12 

4 

+ 3 

4-12 

36 

625 

1 

+ 6 

4 6 

36 

r = 

r 33 











.8 

X 

M 

c*’ 

X 

13 

1 

. +4 

+ 4 

16 

Sums 

37 

... 

-10 

150 













= o*t;6 

14 

0 

+ 5 

0 

0 

Means ( y 2 and y 2 ) 

-'2 7 

4*05 



15 

2 

4-6 

4-12 

72 

<Tp : 


-*v= 

= n/3'98: 

= 2*00 



3 ums 

37 

... 

-28 

410 








Means and u 2 ) 

-•84 

+11 *08 




\ 




0A : 

= Vva 


aA ° - 37 = 

= 3*22 
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8. This requires little explanation. We divide the periods into 
groups of 50 days as shown in the first column, and the number of 
periods in each group is shown in the final column. For the 
values of A no collection in groups is necessary, successive degrees 
being shown in the top line, and the number of corresponding 
cases in the bottom line. In the body of the table are given, in 
large figures, the number of cases corresponding to a definite 
degree for each group of periods. The small figures will be 
referred to presently. 

9. The thick lines indicate the groups which contain approxi¬ 
mately the mean of each co-ordinate. They are chiefly for con¬ 
venience in working (to keep the numbers small), and if a wrong 
group is taken no real harm is done. 

10. The figures below the table represent the whole working. 
We first form v l and v 2 for each co-ordinate by a process which is 
tolerably clear to inspection. For A write down the (top line for 
convenience only—not used in working; and the) bottom line of 
the table. Follow with a sequence of units, with o at the place of 
the adopted mean. Multiply the two columns; and again multiply 
the last two columns. W'e thus get nv 1 and nv 2) where n is the 
number of observations, and v 1 , v 2 two quantities required. 
Dividing by n , we get v x and v 2 , and finally 

°A-( , ' 2 - ,, l 2 ) i== 3 ' 22 ‘ 

[If greater accuracy is required, we should subtract also- 
“Sheppard’s correction” of ^ from v 2 before taking the square- 
root.] 

11. The same procedure with the columns gives us op = 2*00. 

12. We now form the product S(AP). The number in each, 
box is to be multiplied by the little number shown in the corner of 
the box, which will be seen to be the product xy, where x repre¬ 
sents the number of boxes to the right, and y the number up, from 
the adopted mean. The products in the four main divisions of 
the table are shown diagrammatically, and from them we get 

S'(AP)=i 4 2. 

From this we subtract where n is the number of observa¬ 

tions, and , v-^ have been found ; and thus 

S(AP) = S'(AP) - WjVj 1 = 142 -9 = 133 . 

Then the “correlation” r is 

r- S ( AP ) =o- 5 6. . 

n x o- A x o> 

13. This is quite a respectable value for r—anything over o' 5' 
is worth serious attention. But to complete our information we- 
want the probable error of r, which is given by the formula 

, 2?. e = *6 7 (1 — V s ) l fn, 

which can be taken from the following small table (Table III.). 
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For the present case we have n— 37, r = 0*56, and thus probable 
error = ± *08; 


Table III. 

Probable error of r . 


r 


\ 

•I 

*2 

•3 

*4 

1 

*5 

•6 

*7 

*8 

*9 

H 

10 

•208 

*202 

•193 

*477 

•158 

•136 

•107 

•076 

‘041 

20 

•148 

•143 

•137 

•125 

*112 

•096, 

•076 

•054 

*029 

30 

*121 

•117 

*112 

•102 

•09I 

•079 

•062 

•044 

*024 

40 

•104 

•101 

•096 

00 

00 

0 

■ ‘079 

•068 

‘054 

•038 

*021 

-50 

•a'93 

*090' 

•086 

•079 

•070 

•061 

•048 

•034 

*018 

60 

ir> 

C 0 

•083 

•079 

*072 

•064 

•056^ 

•044 

•031 

•017 

70 

•079 

•076 

•073 . 

•067 

XT) 

O . 

•051 

•040 

•029 

•015 

:8o 

•074 

•072 

. *068 

•063 

•056 

048 

*038 

•027 

•015 

90 

•069 

*067 

•064 

•059 

•053 

•045 

•036 

•025 

•014 

*oo 

*066 

*064 ' 

*061 

•056 

•050 

•043 

*034 

*024 

•013 


14. Coming now to the relation between P and A, the mean 
value of A is got by adding v l to our adopted mean, Le . it is 
q°_o°* 84 = 8°*i6,; and the mean value of P is 325-0*27 x 50 — 
311*5 days, the unit adopted in the table being 50 days, by 
which therefore we multiply v v 

Assuming a linear relation between A and P, there are two 
appropriate forms for it, viz. 


(A- 8°*i6) = r. 


o>\ 50 / 


and 


f P ~ 311 ") = r. 5Z(A - 8°-i6) 

V 5 ° > ,°a 


> 


and those new to correlation methods will notice with surprise 
that these are not convertible. To save repeating explanations 
already given in detail elsewhere, we need only remark that to 
find the most probable value of A when P is given is not the 
aarae as to find the most probable value of P when A is given. 
There is the same hind of difference (the analogy must not be 
taken too seriously) as between finding from the fact that a man 
is 25 the date of his probable death, and finding from the date 
of his death the probable date of his being 25 : the latter case is 
clearly complicated by the possibility that he may never have 
reached 25 at all. 

15. In astronomical examples it will frequently happen that 
one of ,the variables, is better determined than the other, and hence 
one pf .the^o alternative forms is indicated as preferable. Thus’ 
in our case P is well known; most of the periods are known within 
& day or two, some perhaps to fractions of a day; and since our 
unit in Table II. was 50 days, the precision is less than one-tenth 
■of the unit. But A may be several degrees wrong, Lei several 
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whole units. Hence it is natural to determine A from P, and the 
appropriate equation .is 


(A - 8°*i6) = o ’56 x 


S'22(F - 311 
5 ° 


or 


2*00' 

A = 2°*6 + *oi8P , 

P being now measured in days , and not in units of 50 days. The 
most probable values of A corresponding to different periods are 
thusi • 1 

For P= ioo d 20o d 3oo d 4oo d 5 °° d 6oo d ' 

A= 4 0 ’4 6°*2 8°*o 9°*8 ii°' 6 i3°*4 


16. This single example will suffice to show the method of 
finding r in such cases. There are many other points to be 
considered^ e.g. how to deal with cases where the relation is not 
linear, or where one of the quantities is not measured but only 
described. But it seems probable that in astronomy the evalua¬ 
tion of r by the above simple process will be the point, most often 
Occurring. 

17. But before leaving the more general topic of correlation 
for the particular application of it we have in view, it is a matter 
of some interest to return now to the simple example given in § 2 ?> 
and work it out as below. 


Table IY. 

The example given in § 2. 


\ A 
B X 

5 

6 

7 

.8 

9 


6 

1 

2 

i 

| 1 

1 

1 



3 

7 


1 

i 

1 


3 

8 


1 

1 

1 

1 

1 

2 

'3 


1 

2' 

3 

2 

1 

9 


«a= 7 I2 = 3 ' 4 6 


°B = ^6=2-45 


A 

Sum. 

Factor. 

vi 

v, 2 

B 1 

3 um. 

Factor. 

n 

v 2 



5 

I 

- 2 

-2 

4 

6 

3 

- I 

-3 

3 

Product S'(AB) 

6 , 

2 

- I 

- 2 

2 “ 

7 

3 

0 

... 

... 

+ 3 

0 

. 7 - • 

• 3 

0 

... 

... 

8 

. 3 

+ r 

+ 3 ; 

3 

0 

+ 3 

1 -8 

2 

+ I 

+ 2 

. 2 

Sums 

9 

. . * 

0 

6 

= 

6 : 

9 

i 

+ 2 

+ 2 

4 

1 







Sums 

9 

• . • 

O 

12 









r = 


9 x J 12 x J6 


= •08. 
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This is a very low value, and on referring to Table III. we see 
that the probable error is greater than o'? : so that the correlation 
of the quantities A and B may be regarded as unsupported. And 
yet at first sight there is a fair appearance of relationship. This 
shows the advantage of having a numerical test. 

18. We proceed now to discuss the interpretation of this and 
similar results. Andin the first place it is an obvious suggestion 
to inquire whether C is also related to P. In a precisely similar 
manner, which need not be given in detail, we find the correlation 
for C and P to be 

r=*30±*io, 

and for the probable relationship 

(C2 7-2) =-o64(P - 311) 

or C = 73 + -ofi4P, 

where C and P are measured in days. But the value of r and its 
probable error show that the probability of a relationship is here 
more doubtful. This is quite possibly due to the much greater 
uncertainty, of the value's of C, which is more difficult to determine 
than A, for reasons which those who have. discussed variable star 
observations will readily appreciate. 

With this formula the values of C would be: 

P=ioo d 200 d 300 d 40o d 5oo d 6oo d 

C= 14 20 26 33 39 46 

19. Now, having obtained the suggestion of correlation and 
two formulae for A and C, I examined in detail some individual 
cases which had come under notice and which are given below. 
It will be seen that there is good confirmatory evidence, so far as 
it goes, of these relationships. 

j 

S Cassiopeia (Ch. No. 432). 

20. In Mem. R.A.S., lv., p. lxix, the Bousdon observations are 
compared with Chandler’s formula in his 3rd Catalogue: 

6io' 5E+ 50 sin (io° E+ 50°). 

From observations E= 17 to 21 (about) a period of 6307 days is 
derived, differing sensibly from 606*2 given by the formula. 
It is remarked : “To obtain [6307 days] from the formula we must 
alter the periodic term in some way, and especially we must 
increase the coefficient, for with the coefficient 50 the greatest 
interval between two consecutive maxima is the greatest value of 

610*5 + 5o[sin 0 - sin (0 4 - io°)] , 

which is 61 o'5 + 100 sin 5 0 or 619*2 days. But we cannot discuss 
this point without reference to other observations. Chandler used 
maxima in 1843 and 1863-93, and any new formula must thus 
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agree with the old one at these periods. Tentatively we may 
suggest some such modification as 

2401590 + 615E + 75 sin (12° E + 90 0 ) 

Chandler's revision gives 

6io*sE + 37 sin (15° E+S9°). 

The value of the periodic term deduced from the period 610 
days by the formulae found is 

47 sin (i3°*6 E + ?). 

The coefficient given by the formula is thus in accord with 
Chandler's first thought, and lies between his revised value and 
the value suggested very tentatively by the Rousdon observations. 
The mean of the three different suggested values, 75, 50, and 37, is 
54; and the fact that the smallest of these, 37, was used in 
deducing the formula accounts for part of the difference between 
the formula value 47 and this mean. 

It is perhaps worth noting that if the suggestion of C=75 
(made quite independently of the present investigation) is 
adopted instead of Chandler’s 0 = 37, the value of r is raised from 
•30 to *43. 

The argument 13°* 6 given by the formula is between the value 
suggested by the Rousdon observations and Chandler's revision. 

R Urssa Majoris (Ch. No. 3825). 

21. In Mem. R.A.S., lv., p. lxxvi, the Rousdon observations 
are compared with Chandler’s 3rd Catalogue formula 

302* iE 4-15 sin (io° E+ 190°). 

It is shown that the Rousdon observations indicate an error of 
150° in the periodic term : so that the coefficient of E should be 
8° instead of io°. Chandler's revision gives 

302'iE + 11 sin (8° E + 238°). 

Pogson's observations of this star have since been very carefully 
discussed by Miss JBlagg; and the outcome of this entirely 
independent discussion (which will shortly be published in 
Mem. R.A.S. , lviii., in the introduction to Pogson's observations) 
was to indicate a periodic term 

40 sin (8° E + ?). 

The formula gives for 302 days a term 

27 sin (8° E + ?), 

in which the coefficient is midway between Chandler's revision 
and the Pogson indication ; and the argument in good accordance 
with both. 

40 
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T Urm Major is ( Ch.‘No. 4511). 

22. In Menu R.A.S., lv., p. lxxviii, the Eousdon observations 
are compared with Chandler’s 3rd Catalogue formula 

257*2E +20 sin (9 0 E +90*) j • 

with the conclusion that the agreement is good. Chandler’s 
revision gives, no change. . The correlation formulae give for a 
period of 257 days a periodic term 24 sin (7* E + ?)' ' 

S Ursx Majoris (Ch: No. 4557). 

23. On p. lxxx the Eousdon observations are compared with 
Chandler’s 3rd Catalogue formula, 

226‘iE+ 43 sin (5°76 E+i8i w, 5) , 

with the conclusion that “ there is a fairly satisfactory accordance, 
though some correction to the formula would improve it.” 
Chandler’s revision gives 

226*5E + 35 sin (5°*4 E + 194 0 ). 

The correlation formulae give for a period of 226*5 days a periodic 
term 22 sin (6°*5 E + ?), which suggests that Chandler’s diminu¬ 
tion of the coefficient has not been carried far enough, though in 
the right direction. 

S Gygni (Ch. No. 7220). 

24. On p. lxxxix the Eousdon observations are compared with 
Chandler’s formula 

322*8E + 15 sin (12° E + 66°). 

It is remarked that the “ periods at (the Eousdon) epoch would 
agree better if the coefficient of E in the periodic term were 
smaller, say 9 0 instead of 12 0 .” 

In Chandler’s revision the periodic term is replaced by a 
secular term 

323E + o*oi5E 2 . 

The correlation formulae give for a period of 323 days 

28 sin (8° E + ?), 

so that the only suggestion made at the time of discussing the 
Eousdon observations is in the direction of better accordance with 
the formulae. 

# 

R Sagittaz (Ch. No. 7257). 

25. This star might have been included in Table I., but it was 
decided to draw the line as regards “long-period variables” at 100 
days. The elements given in the revision are 

7o*56E + 6*5 sin (2**2 5 + 47 Q ). 
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Extrapolating our formulae, we get C^= 12 and A = 3* *9. But if we 
decide that this star may be regarded as a long-period variable, 
we should include it in our table, and then the values of C and A 
would be found closer to those observed, since the star would have 
great weight. Indeed, the correlation for A and P is raised to 
r = o*64± *07, and the formula becomes 

(A-8‘-i)=roi(^p) 

or A= i 0, 9 + ‘020 P, 

which gives A = 3°*3 for P «= 70. 

The value of r for C and P is, however, not much improved, 
being raised from *30 to *32. 

R Cassiop&im (Ch. No. 8600). 

26. On p. xcii of the Rousdon Memoir, the observations are 
-compared with Chandler’s 3rd Catalogue formula, 

429 ' 5 e + 2 5 sin ( I 5 ° E + o°), 
and it is remarked :— 

“The corrections to period shown by the different columns of 
the Rousdon observations are so consistent that it is difficult to 
believe that the mean result can be so erroneous as the formula 
would make it. If the formula were altered to 

431E + 36 sin (i2° E + 9 0 ) 

[there would be a certain improvement]; but it is of course 
impossible to alter the formula definitively without discussing 
other observations.” 

Chandler’s revision, which appeared after the above words were 
in type, gives 

431 *6E + 32 sin (9 0 E + 6o°) . 

The correlation formulae give for the period 431 days a periodic 
term 

35 sin (io°’7 E+ 0 . 

.showing that the improvements suggested by the Rousdon observa¬ 
tions were in the right direction, though not sufficient in magni¬ 
tude; and that Chandler’s revision is in good accord with the 
formula. 

S Delphini (Ch. No. 7431). 

27. Chandler gives no indication of a periodic term in his 
“revision,” printing the period as 277*5 days. In discussing 
Baxendell’s observations (before the present correlation work had 
been undertaken at all), Miss Blagg found clear indications of a 
periodic term at which a preliminary guess of 9 sin (7^° E+ ?) Was 
.made, the 9 being mere guess-work, but the 7^° being indicated witti 
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fair precision. The correlation formulae give for a period 277*5 days 
A = 7**6, which is close to the value guessed. For the coefficient 
C the formula gives 25 days, which is a good deal larger than 
the guessed 9; hut on calculating the coefficient from the separate 
corrections to epoch the coefficient came out 28*3, which accords 
well with the formula. Although the calculation is of a rough 
kind, it seems worth while giving it here to show the kind of 
accordance that is at present obtainable. 


Table IY. ' 
Baxendell’s S.Delphini . 


E 

Observed 
Correction 
to Epoch. 

Argument 
(7J 0 E+const.) 

Sin (Arg). 

Product. 

I. 

- I2 d 

- 26° 

- -44 

+ 5*3 

II. 

~ 5 

- 184 

- *32 

+ i*6 

III. 

+ 3 

- 11 

-•19 

- o*6 

IY. 

- 5 

- 34 

-•09 

+ 0*5 

Y. 

+ n 

+ 4 

+ •07 

+ o*8 

VI. 

+ 20 

+ 114 

+ *20 

+ 0-4 

VII. 

+ 21 

+ 19 

+ ’33 

+ 6-9 

VIII. 

+ 6 

+ 264 

+ ’45 

+ 27 

IX. 

o 

+ 34 

+ •56 

00 

X. 

+ 23 

+ 4 J 4 

+ •66 

+ 1 $'Z 

XI. 

+ 38 

+ 49 

+ 76 

+ 28-9. 

XXIY. 

+ 5 

+1464 

+ ’55 

+ 2*8 

XXV. 

+ 3 

+ 154 

+ •44 

+ 13 

XXVI. 

+ 1 

+ 1614 

+ •32 

+ 0-3, 

XXVII. 

+ 9 . 

+169 

+ •19 

+ 17 

XXVIII. 

+ 6 

+ 1764 

+ '09 

+ 05 

XXIX. 

-19 

+ 184 

- ’07 

+ i* 3 - 

XXXII. 

- 8 

+ 2064 

•- *45 

+ 3*6' 

XXXIII. 

-27 

+ 214 

-•56 

+ 15*1 

XXXIY. 


+ 2214 

.-•66 

+ 8*6 

XXXV. 

-29 

+ 229 

-76 

+ 22*0' 




Sum 

ii8'9 


28. The first column gives the sequence of periods observed bjr 
Baxendell. There is a wide gap between XI. and XXIY. with no 
observations, and it will be seen from columns 3 and 4 that it was* 
unfortunately just in this gap that the corrections to epoch would 
have been largest, and given us an accurate value. ,Still, the. 
periods when the correction vanishes are fairly well determined, 
which is something. They have been assumed to lie midway 
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between IV. and V., and between XXVIII. and XXIX. This is 
only a rough assumption, but the material does not warrant refine¬ 
ments. The “observed correction to epoch” given in the second 
column had been deduced from the observations by Miss Blagg, 
using a purely numerical process,* before there was any idea of 
applying them in this way. In the next three columns the 
product by sin (arg) is deduced, and is almost uniformly positive. 
Dividing the sum 118*3 by 2 sin 2 (arg), which comes out 4'18, ’we 
get as above mentioned 

C= 118*3/4*18 = 28*3. 


S Serpentis (Ch. No. 5501). 

29. For a different reason it is necessary to call attention to 
the case of S Serpentis, period 369 days. Chandler gives C = 116, 
A = 4°, while our formulae indicate 0 = 30, A = 9°. But it is 
apparently possible to satisfy the observations in a different way, 
which assigns to C and A nearly the values of the formula. This 
investigation is given in a separate note following this paper. 

30. What are now wanted are better determinations of C and 
A for stars with long and short periods. Too much depends 
on S Cassiopeise (610 days) and R Virginis (145 days), in the 
present state of our knowledge of them. “The dog is wagged by 
the tail.” But it is the special value of an investigation of this 
kind that attention is directed to special needs which may guide 
observers in selection. The following stars, for instance, seem 
worth special attention from the northern observers:— 


Short Periods. 
Name. No. 

p 

Long Periods. 
Name. No. 

P 

Z Aquilae 

7260 

I27 d 

S Cassiop. 

43 2 

6io d 

W Cygni 

7754 

131 

V Delphini 

7458 

540 

R Vulpec. 

7560 

r 37 

S Cephei 

7779 

486 

R Virginis 

45 21 

r 45 

W Aquilae 

6900 

480 

S Aquilae 

7242 

147 

U Cygni 

7299 

461 

V Capric. 

7571 

157 

Z Sagittarii 

6923 

45 2 

T Herculis 

6512 

i6 5 

R Leporis 

1771 

436 

R Ceti 

845 

167 

T Draco 

6449 

426 

V Tauri 

1717 

170 

R Cygni 

7045 

426 

R Arietis 

782 

187 

R Hydrte 

4826 

425 


31. Returning now to Table I., iu the seventh column is given 
M - m, the interval between maximum and the preceding minimum 
according to ChandleFs revision. For most variables this is less 
than half the period, so that the quantity 

aP = 2(M - m) - P 

* The process can readily be inferred from the description on pp, 66 to 
68 of Mem. R.A.S., vol. lv,, though this description is adapted to a 
period, i.e. a difference of epochs, instead of to the epoch itself. 
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shown in the last column is negative. The quantity a was used 
in M.N ., lxvii. p. 350, in connection with the classification of 
these variables according to their light curves, and it is interesting 
to see whether the type of light curve is related to the quantities 
C and A we have been discussing. The value of r for C and aP 
comes out *oii, shewing that there is complete independence. In 
considering the type of light curve it was early noticed that it 
seemed to he independent of the period,* and it is not surprising 
therefore that C and A, which seem to depend on the period, 
should he independent of the type of curve. The relation between 
A and aP, and those between a alone and C and A, might be 
actually worked out, but it seems improbable that they will yield 
any result of interest. 

32* Column 6 of Table I,, headed “ Phase,” represents a much 
more speculative investigation, but one which it seemed just 
worth undertaking. It. seemed quite possible that C and A, or at 
any rate A, might be directly proportional to P. The limits of 
possible error are such that there is nothing extravagant in this 
supposition, which would be represented by the equation. 

, 7 8*2 X SO 

where fc = -^-==1*32, 

so that A =0264?. 


We may compare the resulting values of A with those of § 15. 


P = IOO d 

200 d 

3 °° d 

4oo d 

s°° d 

6oo d 

A by old equation = 4*4 

6*2 

8°o 

0 

9*8 

O 

n‘6 

0 

T 3*4 

A by new"equation ^ 2*6 

5*2 

7-8 

10*4 

13-0 

15*6 


33. Now what makes this hypothesis of special interest is the 
fact that the cycle in which the term C sin (A° E + const.) is 
completed is 360/A periods or 360P/A days, and if P/A is constant 
we are confronted with the idea of a universal cycle controlling 
all these stars ! The notion seems impossible, and yet there might 
be a simple explanation of it—the cause might be rooted in us, 
and not in the stars. It is difficult to see how any feature of the 
Earth or of the solar system could affect differential comparisons 
of distant stars, but though difficult it may not be impossible. 
The period indicated is 360/*0264 days or 37*4 years. If there 
is a phenomenon of this kind originating with the observer, it 
might have the same phase for all stars, or at least the phase 
might depend on the star’s position. Hence the phases given by 
Chandler were reduced to the common epoch 240 5000, and are 
shown in the sixth column of Table. I. They certainly tend to 
no one value more than another, and a few experiments on their 
being related to a star’s position gave no encouragement. 

34. Another point may. be mentioned. Whether there is 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjournals.org/ at University of Queensland on June 15, 2015 






1908MNRAS..68..544T 


559 


June 1908. of Prof Karl Pearsons Calculation , etc. 


actually a single periodicity or not, the cycle for the periodic term 
is about of the same order for different stars. The value of 
360P/A in years given by the equation of § 15 is 

P=ioo d 2oo d 3oo d 4oo d 5oo d 6oo d 

Cycle = 2 2 y 32 y 37 y 4o y 43^ 46 y 

Now if its coefficient also remains of the same order of magnitude, 
the maximum gradient of the disturbing cause (whatever it may be) 
may be nearly constant. Suppose it constant for simplicity. Then 
the displacement of maximum of a variable is due to the. addition 
to the term 

V\ —/ COS 27 t£/P ' : 

(which represents the ordinary variation near maximum) of a 
term 

2/2= ld . 

where k may be nearly the same for all stars. The time of 
maximum will now be when 


Q= d^ + dy. 


o 27 T j. . 27 rt , 7 

dt ' dp-p f - sm -F +k ' 


or, since t is small, 


h P 2 

4 ?r 2 ./ 


Now / represents the range of variation, and is not very different 
for the stars in Table I., being usually about 5 magnitudes. For 
some stars the range is as low as 3, and for others it is as great as 
7 ; and in a complete investigation the influence of this range must 
undoubtedly be considered. But for the present we will consider 
/ constant. It follows that C, which is the maximum alteration 
of epoch, would be proportional to P 2 rather than to P. Hence it 
is worth inquiring whether the correlation of JC with P is better 
than that of C with P. On working it out, the result came 
r = o*34, which is indeed greater than the r = o*3o previously 
found, but not much. It would seem that for the present our 
best work can be done in getting improved values of C and A for 
individual stars, especially those above given, both by making 
new observations and by collecting and carefully discussing those 
already made. 


SUMMARY. 

§§ 1- 3. Introductory. 

§§ 4- 5. Tabulation of Chandler’s periodic “inequalities.” 

g§ 6—17. An example of working out the “correlation” r 
between two quantities A and P. The result found is r=cr $6 
+ 0*08, which is worth serious attention. The inference is that 
the argument A° in Chandler’s periodic inequalities can be deduced 
from the period P days by the following probable formula:— 

A = 8°*2 + *oi8P . 


§13 includes a table (III.) for the probable error of r. 
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§18. For C (the coefficient of the periodic inequalities) and 
P we find r = 0*30 + 0*10, which is not so conclusive. The relation 
would be 

C= 7*3+ *064?. 

§§ 19-29. Examination of 8 stars in detail, for which special in¬ 
formation was available. The result is favourable. Remarks on 
S Serpentis, which seems to be exceptional. 

§ 30* List of stars needing special attention. 

§ 31. The type of light curve seems to be independent of 
C and A. 

§§ 3 2 ~3 4 * Possibility of a single periodic cause affecting all 
stars considered, but not supported. 


Note on the Period of S Serpentis. By H. H. Turner, D.Sc., 
F.R.S., Savilian Professor. 

1. The forniula given for the maxima of S Serpentis (JSTo. 5501) 
by Chandler in his Revision of the 3rd Catalogue (A.J., No. 553) is 

2388724 + 368*5 E+ 116 sin (4 0 E + 62°). 

The periodic term in this formula attracted attention by consider¬ 
able divergence from the value suggested by the preceding paper. 
The coefficient 116 days is much too large, and the argument 4 0 E 
is too small. With a view to seeing whether it was well established, 
or whether perhaps some other formula would fit the observations, 
inquiry was made of Professor Muller, of Potsdam, who very kindly 
sent a complete list of observed maxima, with full references; 
adding a comparison with Chandler’s formula, the remark that it 
did not fit more recent observations, and a suggestion of his own 
for improving it, which modified both the coefficient 116 and the 
argument 4 0 E in the right direction. But this suggestion does 
not fit the observation of Lalande, and reasons will be given below 
why it is probably too early to suggest a completely satisfactory 
formula, in spite of the fact that Lalande’s observation was made 
in 1794. Hence it seems unnecessary to reproduce here the full 
details, which will doubtless appear in the great work of reference 
for variable stars now being prepared by the German Committee 
of which Dr. Muller is a member. 

2. For our present purpose the observations are sufficiently 
represented by the dates for every fifth maximum shown in the 
second column of Table I. The numeration in the first column is 
that of Chandler. In the third column are given the intervals, and 
we have to decide how to interpret Lalande’s observation, which 
may belong to any one of the epochs - 35, - 34, or - 33. Chandler 
takes - 33; and this gives an average interval of 1876 days for 5 
periods, extending over 7x5 periods iu all. This was quite a 
possible interpretation before the modern observations (represented 
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